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1 Introduction 



An important mathematical problem in Elementary Particle Physics is the evaluation of Feynman inte- 
grals, which usually appear in the perturbative treatment of quantum field theory amplitudes. Besides the 
intrinsic difficulty in solving the integrals associated to a specific graph, in general the number of diagrams 
grows rapidly when the number of loops is increased, which makes necessary to develop methods that allow 
for the automatization both in the generation and the evaluation of such integrals. The first problem that we 
face in dealing with a Feynman diagram is to decide which integral representation is the most convenient in 
order to start the process for finding a solution. Among the different alternatives wc have the parametric rep- 
resentations, in particular the Feynman parametrization {FP) ([1],[2],[3]) and the Schwinger parametrization 
(a-parameters) ([1],[3], [4]) , which allow for the transformation of the loop integrals into scalar multidimen- 
sional integrals. These representations also permit, using dimensional regularization, for a clear and direct 
analysis of the convergence problem, and furthermore the property of Lorentz invariance is also explicit in 
these representations. Recently very efficient analytical and numerical methods for evaluating loop integrals 
have been proposed, which use as a starting point a scalar representation. In particular the Mellin- Barnes [5] 
representation allows for analytical solutions of complicated diagrams, starting from a Feynman parameters 
integral. In numerical calculations an excellent technique is the so called Sector Decomposition ([5], [6], [7]), 
which allows to find the Laurent series of the diagram in terms of the dimensional regulator (e), system- 
atically separating by integration sectors the divergences in the Feynman parameters integral. From this 
point of view it would be convenient to find an accessible way for obtaining the above mentioned parametric 
representations. Although at present there are in the literature algorithms of topological nature [8], its 
implementation is quite complicated from the point of view of the automatization. 

For simplicity we will consider here a scalar theory, although the generalization to fermionic theories is 
straightforward. The final result is a simple algorithm which allows to find the parametric representation of 
any loop integral, and which can be easily programmed computationally. The basis of this formalism is a 
generalization of the completion of squares procedure used in mathematics structures denominated Quadratic 
Forms ([9], [10], [11]) , which are precisely those that appear when applying a scalar parametrization to the 
momentum space integral representations. In essence the expression of the form Q*MQ is a quadratic form, 
where Q is an {L + E) — vector that contains all the independent internal and external momenta of the 
graph and M is a matrix denominated Initial Parameters Matrix (IPM). The end result of the process is a 
recurrence equation that is the support of the algorithm. 

This study is developed making emphasis on the differences that exist in the way of finding the parametric 
representation and the resultant mathematical structure, between the usual method and the one proposed 
here. We also find that the parametric representation can be expressed in two equivalent and directly related 
ways, the first one in terms of matrix elements generated by recursion starting from the IPM and the 
second in terms of determinants of submatrices of the IPM. The relationship between both representations 
is demonstrated in appendices A and B. Finally, two detailed examples are presented, illustrating the 
procedure for obtaining the parametric representation of a Feynman diagram, and which allow to compare 
in practical terms the usual method and the one proposed here. We also add the explicit code to generate 
the recursive elements of the scalar representation, in the symbolic calculation package Maple. 



2 The formalism 

Let us consider a generic topology G that represents a Feynman diagram in a scalar theory, and suppose 
that this graph is composed of N propagators or internal lines, L loops (associated to independent internal 
momenta q = {(^i, g^}, and E independent external momenta p = {pi, ...,pe}. Each propagator or 
internal line is characterized by an arbitrary and in general different mass, m = {mi, ...,mjv}. 

Using the prescription of dimensional regularization we can write the momentum space integral expression 
that represents the diagram in Z) = 4 — 2e dimensions as: 

r = r( \^ [ ^^gL 1 1 

u(p,m) J ^^D/2-i^D/2(^Bf-ml+i5y^-{B%-m% + i6y^- 
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In this expression the symbol Bj represents the momentum of the j propagator or internal line, which in 
general depends on a linear combination of external {p} and internal {q} momenta: Bj = Bj{q,p). 

We also define the set = {ui, un} as the set of powers of the propagators, which in general can take 
arbitrary values. 

Here we will study two well-known parametric representations: the Feynman parametrization and the 
Schwinger parametrization. In the next sections we will show how to express equation (1) in terms of these 
two scalar representations. The technique consists in transforming the product of denominators in (1) into 
a sum through the use of an integral identity. 



2.1 Momentum representation and his scalcir peirametrization 
2.1.1 Feynman Parametrization 

Using the identity: 



N 



N 



— ' ax\...dxN o(l — 2^ Xj)' 



N 



-1 iyi+... + 1'JV 



(2) 



and after defining Aj = (B? — rrtj) , we can replace (2) into equation (1) and thus obtain the following 
generic result: 



G 



n d^'q, 



r(i/i)...r(j/jv) 



J i=i J (i-TT^ 



N N 

j=i j=i 



(3) 



JV 



For simplicity, from now on we use the following notation: die = dx\...dxN H ^ y = {v\ + ... + vj^). 



2.1.2 Schwinger Parametrization 

The fundamental identity for this specific parametrization is given by the equation: 

oo 

-77-= / dxjx]'~'^ e^^{-XjAj), (4) 







which allows, after replacing Aj = (^Bj — mj), to express equation (1) in the following general form: 



L 



oo , \ n dP / \ 

G = ^ / exp E -.-.^ j / exp [- g j ■ (5) 



r(^i) 



The next step is integrating (3) and (5) with respect to the internal momenta, obtaining in this way the 
corresponding scalar parametrization. 
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2.2 Loop momenta integration and pcirametric representation (Usual method) 



The usual way to integrate over internal momenta consists in expanding the sum J2jLi ^j^j reorder 
it in the following manner: 



j=i ' i=i j=i 
or expressed more compactly in matrix form 



N L L L 



N 



Ea;jS/ = q*Aq-2k*q+ J 



(6) 



(7) 



where the following quantities have been defined: 

A =4« Symmetric matrix of dimension L x L, whose elements are functions of the parameters x only: 
A = A{x). 

q L — vector, whose components are the loop or internal 4-vector momenta: q = [qi ...^l]*- 
k =^ L — vector, whose components are linear combinations of external momenta, with coefficients 
that are functions of the parameters x only, so k = k (^x, p) . 

J Scalar term, which is a linear combination of scalar products of external momenta, with coeffi- 

cients that depend on the parameters x only, J = J (x_,p). 

Evidently the specific form of each of these quantities depends on the topology of the corresponding dia- 
gram, and is made explicit once the parametrization formula is applied to equation (1). With the reordering 
presented in (7), we can write both parametrizations and their respective solutions after performing the 
momenta integrations. 



2.2.1 Feynman parametrization 

r(iv.) 



G 



J .7 = 1 J ii-K' 



n d^<i3 



(i7r^/2) 



L r 



q*Aq - 2k*q -|- J - E Xjmj 



(8) 



which once the loop momenta integral are performed, gives finally the Feynman parametric representation: 



N 



dx 5{1 — ^ Xj)- 

3=1 



[detA] 



iV„-(L+l)f 



N 



det A I xjm] -J + k*A-ik 



(9) 



2.2.2 Schwinger parametrization 



G 



1 



/ 



N 



n d^q, 



d X exp ^ XjTnA — J 



\ . 11 " 

/) — 



■ exp (-q*Aq+2k*q) . 



(10) 



In an analogous way, after integration over internal momenta we obtain Schwinger's parametrization of G: 



G = 



y d'x [det A]"* exp 1^ X] ^i''^] - J + k*A-^kj 



(11) 
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The techniques for solving the momenta integrals in (8) and (10) can be found in detail in the literature, 
both for the Feynman parametrization case([l],[2]), as well as for the Schwinger[l] case. This last one is 
usually solved using products of £)-dimcnsional gaussian integrals, in Minskowski or Euclidian spaces. 

Notice that in both parametrizations (equations (9) and (11)) it is necessary to evaluate a matrix product 
that involves an inverse matrix calculation, which makes the procedure not so straightforward to implement. 



2.3 Alternative procedure for obtaining the Pctrametric Representation ( I ) 



Starting from equations (3) and (5), we can choose to represent the term '^jLi XjBj as a function of 

both internal and external scalar products, related through the symmetric matrix M^-'^^, which we will call 
Initial Parameters Matrix (MPI). The dimension of this matrix is therefore {L -\- E) x {L + E). 
For convenience, let us define the momentum: 



Qj 



Pj-L 



if 

if 



L>j>l 
E>j>L 



(12) 



and the {L + E) - vector Q = [Qi Q2 ■■■Q{l+e)Y- 



Using this definition we can reorder the sum J27=i ^j^j^ ^^'^ rewrite it as: 



N 



L+E L+E 



=1 j=i 



(13) 



where M^^^ is clearly a matrix that only depends on parameters. 

The difference in the matrix structure, with respect to the usual method of finding the parametric repre- 
sentation, is that here wc include both external and internal momenta in the same quadratic representation, 
and not only the internal ones as in the usual case wc presented above(see cquation(7)), which produces 
matrix A. In fact, matrix A is a submatrix of M^^^ which already shows an important difference in the 
parametrization starting point, with respect to the usual method. More explicitly we have that: 



/an ... Oil \ 



(14) 



whereas the Initial Parameters Matrix is given by: 

/ an • • • aiL 



1(I, + E) 



a LI 



ClLL 



M, 



(1) 



con = 



(L + E)(L + E) 

if L > z > 1,L > j > 1 
in other cases. 



(15) 



(16) 



In appendix A we show a generalization of the square completion method for diagonalizing Quadratic Forms, 
which is what appears when we parametrize the loop integrals. Looking at the definition of Q in (12), and 
since we have to integrate only the first L momenta, only the main Lx L submatrix has to be diagonalized; 
that is, we need to perform a change of variables in the first L momenta of the {L + E) — vector Q. This 
can be summarized in the following expression: 
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L , _ L+E L+E _ 
3=1 i=L+lj=L+l 

Using the definition (12), the double sum can be expressed in terms of the external momenta: 

L+E L+E E E 

E E Q^M(/^+'^g, = EEM(^^S^+,)^'i-ft- 

i=L+lj=L+l i=lj=l 

Thus the quadratic form of the momenta Q can be written as: 

L , _ EE 



(17) 



(18) 



(19) 



When the square completion procedure is performed to the quadratic form (13), the hnear transformation 
for each internal momentum is given in general by an expression of the form: 



Qj =Qj + f {x,Qj+i,...,QL+E) with j ^ 1, ...,L 



(20) 



(k) 

whose Jacobian is equal to unity. The matrix elements of the type M-j are defined through the following 
recursion relation (see appendix ^4): 



r 



M, 



(fe) 



(fc) 



, if z < (A: + 1) V j <{k + l) 
, in other cases. 



(21) 



kk 



Therefore in a generic way the first L momenta of the vector Q have been diagonalized, using the square 
completion method. Once this has been done, we are in a position to obtain the desired parametric repre- 
sentation. 



2.3.1 Feynman Parametrization in terms of the matrix elements M^j^ 

Using equation (13), the identity (3) can be written in terms of the vector Q, and thus we get the 
following equality: 



G 



r(iv.) 



r(i.i)...r(j/jv) 



J .7=1 J (in^ 



n d'^Qj 



1 



(i7r^/2) 



L r 



N 



Q*M(i)Q - E Xjm 

i=i 



(22) 



Then we expand the denominator of the previous equation, using equality (19), and therefore we obtain a 
more explicit expression with respect to the integration variables Qj: 



G = 



N 



dx 5{1 — E Xj) 



n d^Qj 



(i7r^/2)' 



L 

y 



(23) 



where it has been defined 



N E 



(24) 
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If we now make a second change of variables, such that: 



M. 



1/2 



Qj — > Qj 



M, 



-1/2 



then: 



M. 



io) 



-1/2 



M. 



-D/2 



d Q, con j = 1, ...,L 



(25) 
(26) 



and replacing this in equation (23), we will have the following transformed loop momenta integral: 

L ~ 



G = 



r(i/i)...r(i/jv) 



J .7 = 1 I ITT 



n d^Q, 



M«...Mi^ 



(m^/2)' 



(27) 



In order to perform the integral with respect to the variables Q^- , let us define now the Hiper momentum R 
of (LD) components in Minkowski space, such that: 



d^Q^.-.d^Q^ = d'^'^R. 



LD ] 



Then the expression (27) is reduced to: 

r(A^.) 



G = 



T{vi)...T{vm) 



r ^ f d^^R 

/d^<5(l-E^'.) / — — 1 



(i?2 - A) 



(28) 
(29) 

(30) 



The solution of this integral, with respect to the Hipermomentum R, can be found using the following 
identity: 



- C-lV 



(31) 



and which finally applied to equation (30) gives us the scalar integral, that is the Feynman parametric 
representation of G: 



(-i)^-r(jv.-^) 



G 



N 



d X (5(1 — J2 ^j)' 



j=i 



,2 _ + 



j=i »>i=i 



(32) 



where the matrix elements M^'^'^^^^^^_^.~^ can be easily obtained from the MPI using the recursion formula: 



M, 



M, 



(L+i)L L(L+j) 



{L+i){L+j) - -"(L+»)(L+j) 



M 



(L) 
LL 



(33) 
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(k) 

2.3.2 Schwinger Parametrization in terms of the matrix elements My 

Analogously, using equation (13), the identity (5) can be written in terms of the vector Q as: 

L 

G = ^, , ^, / exp E -^j"'? / r exp - E Q^M^^^Q , (34) 

r(i/i)...r(j/^) y "^y^, ^jj {iT,D/2f f^i ^j' 

or equivalently, using the expansion given in equation (19), we get: 
where again we have defined: 



AT E 

A = E ^jni] - E Ml^tSiL+j) Pi-Pi- (36) 



Now we can solve the momentum integral: 



I f^»p(-|,«fi'0?) - 1 0f«p(-Mffe5).../ 0|exp(-Mffiei). (37, 

In order to find a solution of this integral we make use of the Minkowski space identity: 

J I ITT 2 \ Ji t(7) 2 



which will allow to evaluate (37). Replacing afterwards this result in (35), we obtain finally the Schwinger 
parametric representation for the generic graph G: 

\ ' / 

in terms again of the matrix elements which as we have said before can be readily obtained 

from the MPI using the recursion formula equation (33). 

2.4 Alternative procedure for obtaining the Pctrametric Representation ( II ) 

There exists a direct relation between the matrix elements M>- ' and the determinants of submatrices 
of the Initial Parameters Matrix {MPI). Such a relation can be expressed by the identity: 

M.(''+^^ = -i^ (40) 

'-^kk 

where A-*^^^-' is a determinant which in general is defined by the equation: 
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M['^ ... M« My 



Mki 
Ma 



M, 
Mr, 



kk 



M, 



kj 



M,. 



(41) 



This result is shown in appendix B. There we also present several relations that are fulfilled by these 
determinants and the matrices M^*^^ , and furthermore show how it is possible to evaluate them directly in 
terms of the matrix element Mj^j \ Meanwhile, lot us express the results we have found in (32) and (39), in 
terms of determinants, using for such purpose the identity (40). Then, by direct replacement, we find the 
following final expressions for both parametrizations: 



2.4.1 Feynman Pgirametrization 



After replacing the matrix elements M-'!'^ for the result given defined in (40), the equation (32) is written 



as: 



G 



/N 





^(1) ^(2) 
^00 ■^11 ^{L-l)(L-l) 



N 

E 



E A 
E 



(i+1) 



Pi-Pj 



After a little algebra we get the final Feynman parametric representation: 



(42) 



G 



(-1)^- T{N, - if) 



/N 
(Tx 5{1 — E ^i) 
.7=1 



N E 

^iiEa^j-^i- E ^fLS)(L+j)P^-P3 
j=l »,j=l 



(43) 



2.4.2 Schwinger Parametrization 

In an analogous way, applying identity (40) to (39), we obtain: 



G = 



(-1) 



or simply: 



d X 



A 



LL 



^(1) ^(2) 

^11 ^22 

a(0) .(1)- .(L-1) 
^00 ^11 '^(L_i)(L_i) 



E A^^ + V,,^, 



exp J2 Xj'm'i - E 



A 



(L) 
LL 



Pi-Pj 



G 



r(i.i)..T(i/^) 



dx 



^LL 



exp 



' ^LL E Xj^j - E ^\L+i){L+j) Pi-Pi 

j=l i,j=l 



A 



V 



LL 



I 



(44) 



(45) 



which corresponds to Schwinger's parametric representation. In appendix B it is shown that these determi- 
nants can be evaluated from the matrix elements obtained using a recursion relation starting from the MP/, 
using the following rule: 



(46) 
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This identity is important since it allows to evaluate the determinants that appear in the parametric rep- 
resentations obtained in (43) y (45). We should point out that the matrix A, defined in section 2.2, has a 
determinant which is equal to A^^\ that is: 

detA = Ai^2, (47) 

a very useful identity for comparing more rigorously the different methods for finding parametric represen- 
tations of Feynman diagrams. 



3 The computational code 

The fundamental equation, which allows to evaluate the matrices M^'') starting from the Initial Param- 
eters Matrix is given by the recursive relation: 

<^''=<'--^'tM^^ (48) 
^kk 

o equivalently 

^(fc-i) ^(fe-i) 

^^\k-l){k-l) 

which can be easily programmed in any computer language, and also in a CAS {Computer Algebra System). 
The codification of this equation gives rise to a simple recursive procedure, which we present here in Maple: 

>R:=proc(m,k,i,j) local val: 
>if k=l then val:=m[i,j]: 

>else val:= simplify( R(m,k-l,i,j)- R(m,k-l,i,k-l)*R(m,k-l,k-l,j) / R(m,k-l,k-l,k-l) ): 
>fi:end: 

In this procedure wc have codified the recursive function R(m,k,i,j), where the input parameters are 
given by the following definitions: 

m =^ Corresponds to the Initial Parameters Matrix (MPI), which is obtained at the beginning of the 
parametrization process. Is the matrix that relates internal and external momenta in the quadratic form 
Q*M(i)Q (Q = [qi...qLPi...pE] and m = M(i)). 

k Corresponds to the order of recursion of the matrix. The case k = 1 represents the MPI, and the 
cases k > 1 correspond to matrices obtained by recursion starting from the MPI. 

i,j Is the matrix element to be evaluated. 

The algorithm is very simple. It is only necessary to parametrize the loop integral and recognize the 



r(fc) 

execute the following command or instruction in Maple: 



matrix M^^^ Then we make m = M^^^. Finally if we want to evaluate any matrix elementM|^. , we just 



>R(m,k,i,j); 



4 Applying the algorithm, simple examples 
4.1 Example I 

Now we will compare in actual calculations the usual form and the one presented here for finding the 
parametric representation in terms of Feynman parameters. For that purpose let us consider the following 
diagram, where the masses associated at each propagator are taken as different. 
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Figure 1: Sunset diagram 



First we write the momentum representation of the graph: 



d^gi d^q2 1 1 1 

/ - ml) {Bl - ml) [E'l - m§) ' ^ ' 



where the branch momenta Bj are in this case defined as: 

Bi= qi 

B2= qi+ 92 (51) 
S3 = Pi + q2- 

Applying Feynman parametrization we obtain the foUowing integral: 

1 D D 

G = r(3) J dxidx2dx3 6{l - xi - X2 - X3) J j^^^^—^-^, (52) 


where we define 

n=J2 ^jB] - J2 ^1^]. (53) 

Then, expanding the previous sum and factorizing the result in terms of internal momenta, we get a quadratic 
form in these momenta, which reads: 

3 

= (xi + a;2)(7i + 2a;2 91-92 + (2:2 + 2^3) 92 + 2x3^1.92 + xzv\ - J2 Xjm^.. (54) 



4.1.1 Usual method of finding the parametric representation 

According to the previous formulation (see equation (7)), we can identify the necessary basic elements 
for finding the parametric representation. These are: 

Xi + X2 X2 
X2 X2+ X3 

V ( n V (55) 

k = ( -X3P1 ) 

J = Xspj. 

We start from the general result that we found in equation (9) for Feynman's parametrization: 
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G = 



-i)^-r(iv. - ^) 



/N 
die (5(1 — X] Xj) 
.7=1 



det A Xjm] - J + k*A-ik 



(56) 



In the present case this gives: 



/3 
dxi...dxi5{l- J2 Xj] 



Evaluating the terms that are involved here we get: 



[det A] 



3-3-! 



AT 



det A I ^ Xjm'^j - J + k*A-ik 



3-D • 



(57) 



det A =xiX2 + X1X3 + X2X3 

j^-l _ 1 [ X1+X2 -X2 



k*A-ik 



X1X2 + X1X3 + X2X3 \ —X2 X2 + X3 
xl{xi +X2)pi 



(58) 



X1X2 + XiXs + X2XS 



det A (- J + k* A-^k) = - {X1X2X3) pi 
and considering also the fact that £) = 4 — 2e one finally obtains the Feynman parametric representation: 



1 

G = -r(-l + 2e) Jd-^ 5{1 - EjLi xj)- 



\X1X2 + X1X3 + X2X3\ 



-3+3e 



{xiX2 + X1X3 + X2X3) J2 Xj^] - {X1X2X3) pI 



(59) 



with die = dxidx2dx3. 



4.1.2 Obtaining the scalar representation by recursion 

Remembering the general formula that is used in this method for the parametric representation: 



(-1)^- r(jv. - if) 
r(i/i)...r(i/w) 



G 



d'x 5{l-Yf^^^Xj) 



N ^ (L 1) 

HXjm]- E Pi-Pi 



which in the present case gets reduced to the following: 



G = -r(-l + 2e) / dxidx2dxs 6{l - J2 ^j) 



-2+e 



i=l 



E^jmj-Mg^p^ 



-l+2e • 



From equation (54) one can find immediately the Initial Parameters Matrix (IPM): 

Xl + X2 X2 

M^^) = M = I X2 X2 + X3 X3 
X3 X3 



(60) 



(61) 



(62) 



12 



It is now only necessary to calculate the matrix elements using the recursive function described in section 
3. Basically we need to evaluate the following identities: 



= i?(M, 1,1,1) 
M^^ = E(M,2,2,2) 
i?(M,3,3,3). 



M. 



(3) 



33 



The results, after writing the commands in Maple, are respectively: 
>R(M,1,1,1) ; 



>R(M, 2,2,2) ; 



>R(M,3,3,3) ; 



Xi + X2 



X\ + Xi 



X1X2XS 



XIX2 + XIX3 + 2:2X3 

Thus replacing these expressions into equation (61), we obtain: 



G = -r(-l + 2e) j d-i (5(1 - Xj\ 



\X1X2 + X1X2, + X2Xs\ 



-2+e 



E 



XiX2X:i 



1 i 

j = l ■' X1X2 + XiXs + X2X3 



-l+2e ' 



and then we have the same scalar representation as found before in (59): 



(63) 



(64) 



G = -r(-l + 2e) J dl^ S{1 - EjLi Xj) 



[X1X2 + X1X3 + X2X3\ 



-3+3e 



H -l+2e • 



{xiX2 + X1X3 + X2X3) J2 Xjrn] - xiX2X3p'^ 



(65) 



4.2 Example II 

Let us consider now the following diagram: 
The loop integral is given in this case by: 



G 



iTrD/2 [Bl - ml) (B| - m2) (B| - ml) ' 



where the branch momenta Bj have been defined in the following way: 

-Bi =(?i 

-B2 = pi + gi 

B3 =Pi +P2 + qi- 

The next step is to apply Feynman's parametrization, obtaining the following integral: 

d^gi 1 



G = r(3) j dxidx2dx3 5{1 - EjLi Xj) j 



(66) 



(67) 



(68) 
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Figure 2: Triangle diagram 



where the denominator Q is given in terms of the internal momenta by: 

fl = {xi+X2+ X3)ql + 2 [{x2 + xz)pi + X3P2] -qi + {x2 + 2:3) + 2x3Pi.p2+ 
X3Pl~Tf']=iXjm'^j. 

4.2.1 Usual method of finding the parametric representation 



(69) 



Starting from equation (69) we can recognize right away the basic necessary elements for finding the 
parametric representation. These are: 



A = (a;i + 2:2 + 2:^3) 

k = -(x2 + 2:3)^1 - X3P2 

J = [X2 +X'i)p\ + 2X3PI.P2 +X3pI, 

and therefore the resulting scalar integral will be given in this case by the expression: 



G = (-l)3r(3 - ^) / dxi...dx3 5(1 - i: X,). 



[det A] 



3-D 



N 



det A ( ^ Xj-m^ - J + k*A-ik 



Evaluating each term, we obtain: 



det A 

A-i 

k*A-ik 



=xi + a;2 + ^3 
1 

Xi + X2 + X3 
_ [{X2 + X3)pi + 3:3^2] 
Xi + X2 + X3 



and 



detA(-J + k*A-ik) 



-{xiX2 + XiX3)pI - 2xiX3pi.p2 - {X1X3 + X2X3) pI 

-XiX2p\ - X2X3PI - XiX3{pi + P2f 



-XiX2p\ - X2X3PI - XiX3pI, 



(70) 



(71) 



(72) 



(73) 



where we have used the condition {pi + P2)'^ = P3, and then put D = 4 — 2e. Thus we finally arrive at 
Feynman's parametric representation: 
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G = -r(l + e) I dx (5(1 - E 2;^)- 



[Xi +X2+ Xz] 



-l+2e 



+ a;2 + 0:3) X) Xjjv?, - {xiX2)p\ - {x2X3)pl - {xiX3)pl 



, (74) 



withrfa; = dxidx2dx3. 



4.2.2 Obtaining the scalar representation by recursion 

In this method the general formula for the parametric representation is: 



G = 



/ die (5(1 — ^ Xj)- 
J j=i 



which in our case is reduced to the following in D = 4 — 2e dimensions: 



N E 



(75) 



G = -r(l + e) / dxidx2dx3 S{1 - T, Xj)- 



M- 



-2+e 



E x,m2 _ m^^pI - Mg^pi.p2 - Ml^'p2.pi - M^^'pI 



/-(2)„2 



(76) 

The next step consists in the evaluation of the matrix elements of ' . In order to do this, and starting 
from equation (69) we find the Initial Parameters Matrix ( MPI): 



M(i) = M 



Xi+ X2 + X3 X2 + X3 X3 
X2 +X3 X2+ X3 X3 

X3 X3 X3 



(77) 



Using the recursive routine proposed in section 3, the necessary matrix elements M^^ are evaluated: 

m[1^ = R{M, 1, 1, 1) 
M^^ = R{M, 2,2,2) 
M^f = M^^ = R{M,2,2,3) 
M^^ = R{M,2,2,2), 

and executing the Maple commands, we get the following results: 
>R(M,1,1,1) ; 

X1+X2 + X3 



(78) 



>R(M,2,2,2) ; 



X1X2 + X1X3 
X1+X2+ X3 



>R(M,2,2,3) ; 
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>R(M,2,3,3) ; 



X1+X2+ X3 



X1X3 + X2XS 

Xi + X2 + X3' 



For the sum J2tj=i ^(L+J(z,+j)Pi-Pj, we get 



^ ^AL+l) _ {xiX2+XlX3)pl + 2xiX3Pi.p2 + {xiX3+X2X3)P2 

,2^/^\L+i){L+j)Pi-Pj 



■ X3 

{xiX2)Pi + {x2X3)pI + {X1X3)P3 



Xi+ X2+ X3 



Thus, replacing these quantities in (76), we obtain: 



G = -r 



(1+.)/ 



dxidx2dx3 5{1 — X) Xj 



A 2 _ ixxX2)pl + {X2X3) pI + {XIX3) pI 



X1+X2+ X3 

which finally is reduced to the same parametric representation deduced before in (74): 



G = -r(l + e) I dx (5(1 - Xj) 



[xi +X2+ X3] 



-l+2e 



{xi +X2 + X3) J2 XjUij - {X1X2) pI - {X2X3) pI - {X1X3) pI 



where die = dx\dx2dx3. 
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Summary 



Usual Parametric Representation 

/ an 



aiL 



a-LL J 



Q = [Qi 12 =^ Only internal momenta. 



FEYNMAN PARAMETRIZATION 
Before 



T{i.i)...TMi ^ ''J (i^. 



n d^qj 



N 



q*Aq — 2k*q + J — Yl 



After 




(_ipr(jv^_j^) 

, . , — / dx (5(1 - Xj) 



SCHWINGER PARAMETRIZATION 
Before 




After 



(-1) 
T{ur)...T{yN) 



I n °° / 
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Alternative Parametric Representation Form in terms of the matrix ele- 
ments m/*^^ 



\ (l+e)! (l+e)(l+e) / 

Q = [91 Q2 P1P2 "-Pb] => Internal and external momenta. 



a- 



FEYNMAN PARAMETRIZATION 
Before 



G 



r(iv.) 




After 



(-1)^- T{N, - if ) 

rK)...r(i/jv) 



N 



dx 5{l - J2 



,2 _ 



E^jm]- E Ml^+i)\L+j) Pi-Pi 



SCHWINGER PARAMETRIZATION 
Before 



ft rf^'Q. 



—— — 7- / da; exp > a;,mf / 



N 



■exp - E Q*M(i)Q 



After 



G 



d X 



■° I N E 

,1 ;\/r(-f'+l 



\J=1 »J=1 
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5.3 Alternative Parametric Representation Form in terms of subdeterminants 
of the MPI 



FEYNMAN PARAMETRIZATION 
Before 



G 



nNu) 



T{uy)...T{vN) 



r N 

J 7 = 1 




After 



G = 



i 

/N 
(Tx (5(1 — Xj) 



SCHWINGER PARAMETRIZATION 
Before 




After 



G 



exp 



f ^ 2 A(i'+1) \ 



A 
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6 Conclusions 



There are two main aspects that need to be emphasized in the present work. The first is the simplicity 
of the method, both in the actual calculation and in its application to a particular topology. From the point 
of view of the mathematical structure of the final scalar representation, there is a remarkable difference 
with the usual method. In the usual parametric form of a loop integral, it is necessary to evaluate a scalar 
term and a matrix product that involves an inverse matrix calculation. The method proposed in this work 
is based on a simple change in the initial procedure in the search for a parametric representation of the 
momentum integral, so that both the scalar term and the matrix product with inverse matrix are included 
in an expansion of internal products of external momenta, in which the coefficients of such expansion are 
determinants of submatrices of the matrix that relates internal and external momenta (IPM). The evaluation 
of determinants is much simpler than the evaluation of matrix inverses. Moreover, the most important aspect 
is that such determinants can be in turn calculated from matrix elements obtained using a recursion relation 
starting from the IPM, in a simple and straightforward way. 

The second relevant aspect is that this method can be easily implemented computationally. This allows 
for a fast automatization of Feynman diagram generation, obtaining simply and directly the parametric 
representation as a step towards a complete numerical or analytical evaluation whenever possible. 
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APPENDICES 



A Quadratic Forms and its diagonalization by square completion 

A quadratic form in n variables is an expression which can be written in matrix form as the product 
x*Mx, where x is an n-dimensional vector, given by x = [xi, Xnf, and M is a generic nx n dimensional 
matrix. That is: 

n n 

x*Mx= E Ea^i^ij-a^j- (82) 
i=i j=i 

Let D be an n X n diagonal matrix. The expressions x*Mx and y*Dy are equivalent if there exists a linear 
transformation y = Px and y = Px such that x*Mx — y*Dy, that is: 

M = P*DP. (83) 
The quadratic form is then transformed into a sum of n linear terms of the type y^yjSij. 



A.l Squaire Completion Procedure 
A. 1.1 Completing the square for xi 

Every Quadratic Form can be diagonalized using the square completion procedure, which generates the 

required linear transformation. First wc define a matrix M = M^^\ and then wc expand the matrix product 
x*]V[*^^)x in order to complete the square associated to the parameter xi. Then we arrive at the following 
result: 



E E ^^MlPxJ 

2=1 j=l 



\3=2 



Ij 



E x.M^p ) a;i + E XiM'^'x^ 



i=2 



ij=2 



xl+xi\J2 



M- 



(1) 



M 



E X- 



, i=2 Al- 



ii 



i=2 M^l\ 
Let us define now the new variables: 



XI + E 



n M. 



:i) 



(1) 

11 , 



Xl 



E X^M,, 



i,j=2 



ij -^3 



•-'Mir', 



(1) 



2/1 



yi 



Xl 



E 



M. 



(1) 



Xi 



il 



M. 



(1) 



i=2 ^"11 

Xl+2^ {1)^3^ 
j=2 -'«11 



and also the matrix M^^' = IMj^^''!, such that the elements of this be given by the relation: 



(1)^1^(1) 



M},>M. 



M. 



(1) 
11 



(84) 



(85) 



(86) 



Therefore we can rewrite the quadratic form (82), with the first parameter already diagonalized, in the 
following way; 



i=2 j=2 



ij Xj- 



(87) 
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The second term in the right hand side can be simphfied, since from the equation (86) one obtains that 
Mji^ = M[f = 0, for l<j<n. Thus we can write: 

E E XiMg^Xj = E E XiMg^xj = x*M(2)x. (88) 

i=2j=2 1=1 j=l 

In summary, in the quadratic expansion the first term has been aheady diagonahzed, a fact that can be 
described by the foUowing expression: 



(89) 



A. 1.2 Completing the square for X2 

Now we take the second term of (89), and the same procedure foUowed above is repeated in order to 
complete the square for the parameter X2, which gives: 



X*M(2)x = E E X^M, 



j=2 j=2 



ij -^3 



M^^xl+X2 [ E M^^^xA + (e xMf] + E XiM^fxj 

\i=3 / \i=3 / j,i=3 



Xi+Y.X^ 



(2)^ 



i2 



n M. 



i=3 M. 



(2) 



X2+ E 
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=3 M. 



(2) '"J 
22 / 



E Xi\Ml.' I a;,-. 



i.j=3 



M. 



(2) 
22 



Let us define, analogously to equation (85), the new variables: 



X2 + E 



(2) 



M! 



y2 = X2+Y. 



i=3 '"22 
n ^(2) 



j=3 ^22 



and the matrix M^^) = |m^?^^ |, where we set M^f = M.^J 



(2) 



M. 



(2) 
22 



. Then we obtain: 



(90) 



(91) 



x*M(2)x = Mi^y^ya + E E XiM^fxj. 



The expression for the matrix element m\^^ implies that M^-' = — 0, with 2 < j < n, and since we 



i=3 j=3 



f(3) 



r(3) 

^i2 



(92) 



^3 

also had that M^^^ = Afj? = 0, then M^) = uf^ = 0, where I < j < n. In this way we can write the 
second term as: 

E E XiMlfxj =J:E XiMlfxj = x*M(3)x, (93) 

z=3j=3 i=\ j=\ 



and therefore now the first two components of x have been diagonahzed: 

x*M(i)x = M^iViyi + M^^^y^yi + x*M(3)x. 



(94) 
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A. 1.3 Generalization of the Square Completion Procedure for x. 



Notice that the last term in (94) is another quadratic form, which then will allow us to complete the 
square for the parameter x^. The procedure can be repeated successively for X3,...,Xn, and therefore the 
following relations are determined by induction: 



yi = xi + EHi+i X: 



(0 



(0 



yi = xi+ 2^ . 



M, 



(0 







(95) 



M} 



(0 



M, 



if i < (/ + 1) V j <{l + l) 
in other cases. 



Here the matrix M^'') (con 1 < fc < n) has the following generic structure: 

/O ••• \ 



(fc) 
fefe 



M. 



kn 



\ 







(96) 



In general, the procedure of square completion of the k-element of x, for 1 < fc < n, transforms the initial 
quadratic form into: 

x*M(i)x = M^i'yij/i + ... + Mll^y^yk + x^M^^'+^^x. (97) 

The complete process, that is after n square completions, diagonalizes the quadratic form x*M(^^x and 
transforms it into a diagonal bilineal structure, of the form: 



x*M(i)x = M[l^y,y^ + ... + M^^y^yr, = y*Dy, 



where we identify: 



D = diag 



y = [yi,-,ynf 



— 1* 



(98) 



(99) 



A. 2 Some Properties 



1. The relation between the vectors y and .t, is defined by the equation (95), and from it we can identify 
the transformation matrices that fulfill the equations: 



y = Px A y = Px. 
Specifically it is possible to determine P and P, given by: 

/ 1 



P = 



\ 



(100) 















1 











1 



M, 



(101) 



1 / 
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( 1 



V 

2. From the equations in (95), we find that: 



(1) 



M, 



M("+i) = {0} 

3. The transformation matrices P and P have the following property: 

det P = det P* = 1 
detP = detP* = 1 



4. If M^^^ = [M^^^]* (symmetric case), then the following identities hold: 



p = p 

x*M(i)x = y*Dy 
]V[(i) =ptDP, 



where the matrix D is the diagonal matrix given by: 



D = diag 



M[l\M^f,...,M^l^ 



A. 3 Evaluation of the determinant of M^^^ 

From the previous results, the determinant of M^^^ is given by: 

detM(i) = detP*DP = det P*. det D. det P = M^l^ M^l\..M^l\ 
The conditions for evaluating this determinant are given in Appendix B. 

B Matrices M^'^) 

B. l Generalization of the matrices M*^*^) 

It is possible to generalize the n x n dimensional matrices M^'^-' starting from the recurrence 



(fe) 
fefe 



As an example let us consider a generic matrix A^xn = {O'ij}^ and define an input matrix M^^^ = 



24 



B.1.1 Generating M^^) 

Let us evaluate the particular cases of the first row and first column. That is: 



'11 ^"ij 



0, 



= 0, 



(1) 
11 



(i = l,..,n). 



(109) 



(110) 



The other matrix elements do not present a particular interest, and are evaluated using the recursion relation 
(108). Then the matrix M*^^' gets structured in the following manner: 



/ 



M. 



(2) 



22 



\ 

.(2) 
2n 



V M, 



(2) 
n2 



(111) 



Notice that M^^) is computable only if M]^' ^ 0. 

B.1.2 Generating M^^) 

Having M^^^ already evaluated one can construct M^^^ Let us analyze the first and second row. For 
the first row we have that: 



'ij - ^"ij 
while for the second row: 



Mff=M(?- ^^^r^ =0, 



M. 



(2) 



since M^^ = 0. 



22 



{j = 1, ...,n), 



(112) 



Mi'?=MiV- J' =0. 



^2j 



'2j 



(2) 
22 



Analogously, for the first and second column we have the following values respectively: 



(2) 



M. 



(2) 



0, since M.\ 



(2) 



{i = l,...,n) 



and 



22 



(2) 



(2) 



0. 



22 



(113) 



(114) 



(115) 



The other elements have values according to (108). Finally the matrix M^^^ gets the following form: 

/o X 











M. 



(3) 
33 



M.. 



(3) 



3n 



(116) 



\0 m!^^ ■■■ M^^2. J 



(2) 

The matrix is defined only if the matrix element 7^ 0. The procedure can be repeated successively for 
the rest of the matrices generated by recursion, thus finding that for k G [1, n] one gets: 
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/ 



V 



(fe) 
fefe 



with the condition that M^'^) is defined only if M],^ 7^ or fc = 1, 2, n - 1. 



M, 
(fe) 



fen 



\ 



(k) 

kn 



(117) 



B.2 Elements M}: 



(k) 



From the previous results we van find the relation that exists between the matrix elements generated 
by recursion and the input matrix elements M^^^ = A„xn = {^ij}- For M-f ^ 



M^^W*^) M^^W^^) M^^^M^^'> 



1] 



o equivalently 



For Mff> we have that: 

''J 



Mil' 



an aij 



an 



(2) _ Mg'Mg^ _ 



ail aij 

a^i ciij 



an 



Some simple algebra gives the following result: 



M. 



ail ai2 




an 




an ai2 




021 




|aii| 


ail ai2 
a2i a22 





(118) 



(119) 



(120) 



an ai2 ay 
a2i 022 a2j 

an ai2 aij 



an ai2 
a2i a22 



Let us now define the determinant A-*^"*^^^, such that it corresponds to the determinant of a submatrix of 
the input matrix M^^^ = A„xn, whose dimension is (fc + 1) x (A; + 1), and which is given by the following 
identity: 



(121) 



,(^+1) 



an 

aki 
an 



aik aij 



akk akj 
a%k aij 



(122) 



Let us see the following examples: 
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B.2.1 Example I 



,(2) 
^34 



ail ai4 
034 



(123) 



B.2.2 Example II 



^(3) 
^33 



ail 0-12 ai3 

021 O22 ^23 
0131 «32 033 



(124) 



Applying this definition in the equations (119) and (121), we obtain: 



(2) 



A 



(2) 



A'- 
^11 



M, 



(3) 



A, 



(3) 



,(2)- 



(125) 



(126) 



Through an induction process we can directly generalize the relation that exists between the matrix elements 
of M-j^ and the input matrix M^^^ = A„xn- In general one gets: 



M 



(fe+i) 



A^" 



(fc+i) 



^kk 



(127) 



B.3 The matrix M^*^) in terms of determinants of submatrices of M*^^) 

In appendix A it was previously shown that the determinant of the input matrix M^^^ = A„xn is given 
by the expression: 



det A = detM^i) = M^^ M^^ ...M^^ . 



(128) 



Using equation (127) we can write the matrix elements M^f,' as ratios of determinants of submatrices of 
M(i). Then we have that: 



^11 ^22 



det A„xn = detM^^) 



a(0) a(i) "■" a("-2) a("-i) 

'^OO ^11 ^(n-2)(n-2) ■^(n-l)(Ti-l) 



a(") 

^00 



(129) 



Here Aqq^ = 1, which can be shown by calculating the determinant of a scalar. Let us evaluate the determi- 
nant of the matrix M^^^ = Aixi = (aii), that is: 

det Aixi = detM^^) = an- 



On the other hand we have that: 



det Aixi = detM^i) = m[1> = 



A 



(0) ■ 

00 



(130) 



and applying equation (122) one obtains that A^^,^^ = an, which by comparison gives: 



A, 



(0) 

00 



1. 



(131) 
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Finally it is shown that: 



detA„xn = detMW = Ai") 



an 



air 



(132) 



a result that is evidently correct. In summary we can rewrite the matrix M^*^^ in terms of subdeterminants 
of M(i) = A„xn, that is: 



/ 



,(fe-i) 
\fe-i)(fe-i) 



\ 



A 



(fe) 

kk 



A 



(k) 

kn 



\ 







(133) 



Notice that the relation of the recursive matrix elements with the ratios of determinants provides the con- 
dition for evaluating the matrix M*^'^^. This is that the determinants A||^_Jj^j^_j^^ (Principal Minors) be 
non-vanishing, a condition that is evident in identity (133). 

B.4 Evaluation of determinants A,-'^ in terms of the matrix elements of M^*^^ 

The relation that we found for the recursive matrix elements in terms of a ratio of determinants is given 
by the equation: 



M, 



A(fe+1) 
(fe+1) _ ^ij 

^fefe 



(134) 



We can reorder this such that: 



'^ij — ^kk ^"ij ' 



(135) 



where A[,^' corresponds to a determinant called the Principal Minor of order k x k, which can be expressed 
directly in terms of recursive matrix elements, such that: 



A«=M|i\..M^, 
and therefore we obtain the identity: 

which allows for the possibility of evaluating any subdeterminant of the matrix M^^^ . 



r(i) 



r(fe) 



(136) 
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